Graphene in periodically alternating magnetic field: unusual quantization of the 
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We study the energy spectrum and electronic properties of graphene in a periodic magnetic field 
of zero average with a symmetry of triangular lattice. The periodic field leads to formation of a 
set of minibands separated by gaps, which can be manipulated by external field. The Berry phase, 
related to the motion of electrons in k space, and the corresponding Chern numbers characterizing 
topology of the energy bands are calculated analytically and numerically. In this connection, we 
discuss the anomalous Hall effect in the insulating state, when the Fermi level is located in the 
minigap. The results of calculations show that in the model of gapless Dirac spectrum of graphene 
the anomalous Hall effect can be treated as a sum of fractional quantum numbers, related to the 
nonequivalent Dirac points. 

PACS numbers: 73.21.-b,73.50. Jt,75.47.-m,73.23.Ra 
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I. INTRODUCTION 

The energy structure of graphene includes the Dirac 
points where the energy spectrum of electrons is rela- 
tivistic. This can lead to many very unusual transport 
properties of this material^"— . One of examples is the un- 
usual integer quantum Hall effect in graphene. As found 
theoretically and confirmed experimentally^'^ there is no 

= plateau, which is related to the Dirac point in the 
electronic spectrum. The other very important property 
is that the Landau splitting is very strong in graphene, 
which makes it possible to observe the quantum Hall ef- 
fect at room temperatures^. 

Recently, we proposed to use a periodic magnetic 
field for quantization of the electron energy spectrum of 
two-dimensional electron gas^. Our calculations demon- 
strated that periodically alternating field leads to an en- 
ergy spectrum with a number of minibands separated by 
energy gaps. In such system, the anomalous Hall effect 
(AHE) is nonzero, and it is quantized if the chemical po- 
tential is located within the gap. It was also suggestec^^S, 
that the realization of such structure under periodic mag- 
netic field can be easily achieved by using a lattice of 
magnetic nanorods*^. It should be noted that another 
possible way is to use the 2D skyrmion lattice'^*, which 
have been recently observed in thin layers of helical mag- 
nets (however, in this case the temperature should be 
rather low)^. 

Graphene has obvious advantages to be used instead of 
semiconducting quantum well with two-dimensional elec- 
tron spectrum. Graphene is naturally two-dimensional, 
and the technology of graphene-based structures is much 
simpler. Besides, one can expect much stronger effects 
related to the magnetic quantization due to the specific 



parameters of graphene. 

In this work we calculate the electron energy spectrum 
and the Chern numbers characterizing topological prop- 
erties of the energy bands. The results are essentially 
different for graphene as compare to the 2D electron gas 
with parabolic energy spectrum. We find that the Chern 
numbers of the bands corresponding to higher-in-energy 
excitations of electrons and holes are zero at sufficiently 
large field. On the contrary, the low-energy bands of 
electrons and holes have nonzero Chern numbers, which 
results in the quantized anomalous Hall effect. 



II. GENERALIZED MODEL WITH THE GAP 

We start with a generalized model, which includes the 
relativistic Hamiltoniar*^ describing the electron energy 
spectrum near the Dirac points K,, K,' in the Brillouin 
zone of graphene (Fig. la) and an additional term leading 
to the gap 



H 



K,K' 



ieA, 



-At,, (1) 



where the components of vector r are the Pauli matrices 
acting in the space of graphene sublattices. Introducing 
the gap parameter A makes the A and B sites different 
in the crystal lattice of graphene. Physically, it can be 
the case of different atoms in the honeycomb lattice like 
in two-dimensional boron nitride. The model with A ^ 
can be also realized for graphene in periodic electric 
and magnetic fields discussed in Ref. [ij]. In the case 
of graphene we should put A = but for pedagogical 
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reasons it is instructive to start from A 7^ 0. The vector 
potential A(r) describes the effect of magnetic field. 

We assume that the graphene sheet is in the periodi- 
cally alternating magnetic field in z direction. In our an- 
alytical calculations wo consider one-harmonic approach, 
S(r) = (Bo/3) Rc^''^^ e'S' "', where the in-plane vec- 
tors gj can be viewed as basic vectors of the triangular 
lattice. In the numerical calculations we include more 
harmonics to make the consideration realistic, preserv- 
ing the triangular symmetry of the field-imposed lattice. 
One can use the gauge in which the correspondent vector 
potential is also periodic, A(r) = Im^^ a; e*^* ''^ where 
Bj = (i?o/3(??) (no X gj) and no is the unit vector along 
axis z perpendicular to the plane. 

We also assume that the period a of alternating field is 
much much larger than the lattice constant of graphene 
ao- The large scale field-induced triangular lattice deter- 
mines the corresponding small Brillouin zone in the in- 
verse space (Fig. lb), and we denote the symmetry points 
of this zone by K, M, etc. It should be stressed that these 
points are not related in any way to the symmetry points 
of the Brillouin zone of graphene, which we denote like 
K and K' (Fig. 1). 



III. PERTURBATION THEORY 

Let us asstime A very small and positive, A > 0. In the 
case of weak periodic field, the term with A(r) in (1) is a 
weak pcrtinbation. We consider first the states related to 
the Dirac point /C. For A = 0, the unperturbed solution 
for the low-energy state in the F point with e > is 
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i(r) = 



1 



(2) 



For the lower-in-energy states with e > in M and K 
points we find 



V'M,-M(r) 
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Using these functions we calculate the Berry phase 7c 
along the contour V-M-K-V, shown in Fig. lb, as 



7c = arg(r|M) {M\K) (/^|r) = 0, 



(6) 



Thus, for the contour along the whole Brillouin zone 
and the first energy band we obtain 7+1 = 127c = 0. 
Correspondingly, the Chern number of the (+1) band 
is Ch+i = 7+i/27r = 0. Using the same method for 
the negative energy band (—1), we find 7_i = — 27r and 
Ch_i = -1. 






FIG. 1: (a) Brillouin zone of graphene with two nonequiva- 
lent Dirac points; (b) Brillouin zone corresponding to a single 
Dirac point of graphenc in a large-scale periodic field of tri- 
angular symmetry. Contour 7c in the fe-space. 



where S the sample surface, and index ±1 refers to the 
states in the first low-energy energy band with posi- 
tive and negative energy, respectively. Assuming A very 
small, for any other point k ^ the solution is 
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2S V + iky)/k 



(3) 



corresponding to the states with energy ek.±i — ±?;fc. 

The Hamiltonian of interaction is Hint = ''' ' 
Matrix elements of the periodic field between the states 
^ka ^^"^ V-'k'^ 9,rc nonzero only for k — k' = g, where 
g is one of the vectors ±gi, ±g2, ±g3. Correspond- 
ingly, the weak periodic field mostly affects the states 
at the Brillouin zone edge, which have close in energy 



(0) 



-(0) 



counterparts, e^^ ~ ^k-g 



Then for the state at the 

F point V'r(r) — V'r'''(i"): whereas the wave functions 
near the points M and K of reciprocal lattice can be 
presented as a superposition of unperturbed functions. 



V'k(r) 



coV'k^(r) 



■E,c.<i,,(r). 



Wc should also take into account the contribution of 
the another non- equivalent Dirac point /C'. For this point 
the Hamiltonian (1) differs by the opposite sign before 
dx and A^. Calculating the wavefunctions in this case 
and using (6) we can find for the positive energy band 
Ch_|_i = — 1 and for negative band Ch_i = 0. 

If wc take A < 0, we obtain for the JC point Ch+i = — 1 
and Ch_i = 0. Correspondingly, for the fC' point the 
results are Ch+i = and Ch_i = — 1. 



IV. SYMMETRY ARGUMENTS 

These results can be understood using the symmetry 
of Hamiltonian (1). Indeed, the transition from /C to /C' 

with simultaneous reversion of the sign of A is related to 
the symmetry under the unitary transformation (rotation 
in sublattice space) 



k-gi 



r-'H'^{A)Ty = H'^ (-A) 



(7) 
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Such transformation applied to the wavefunctions does 
not change the Berry phase (6). Therefore, the corre- 
sponding Chern numbers should be equal, 

Ch„(/C,A) = Ch„(/C',-A). (8) 

Besides, changing the sign of field (i.e., A — > —A) in (1) 
with simultaneous spatial inversion and reversion of the 
sign of A is equivalent to the change of sign of Hamilto- 
nian 

•H(r,A,A) = -■H(-r,-A,-A), (9) 

which corresponds to the electron-hole symmetry. It can 
be also understood as a symmetry to the pseudotime in- 
version, T = iTyTZ, where TZ is the complex conjugation 
operator 

f-^HiA,A)f ^H{-A,-A). (10) 

The complex conjugation operator acting on the wave- 
functions changes the sign of the Berry phase, and we 
obtain 

Ch„(A, A) = -Ch„(-A, -A). (11) 

Also, there is a symmetry with respect to Tz transforma- 
tion: 

T-'H{A)rz^-n{-A) (12) 

leading to Ch„(A) = Ch_„(-A). 

Using these properties, we can analyze what happens 
when the parameter A changes sign from negative to pos- 
itive. It will helps us to determine properly the AHE for 
the case of A = 0. 



V. ANOMALOUS HALL EFFECT 

We start from relation between the Chern numbers 
and AHE. For the contribution of points /C,/C' to the 
off-diagonal conductivity, in the case when the chemi- 
cal potential fi is in the gap, one can use the standard 
formula^- 

O'xy — 2^ /(ekn)f^kn, (13) 
kn 

where /(e) is the Fermi-Dirac distribution function, 
f^kn ~ Vk X is the Berry curvature, and = 

— i(krt|Vk|kn) is the gauge connection. The energy spec- 
trum is shown in Fig. 2 for different values of magnetic 
field (the discussion is in Sec. VII). 

The sum over n in Eq. (13) formally includes infinite 
number of energy bands with negative energies of elec- 
trons. However, one can show that in the model of Eq. (1) 
the contribution of all the bands with £kn < to the AHE 
is equal to zero. The reason is that if we take A ^ and 
the chemical potential /i = than the sum of Chern num- 
bers for the /C and /C' Dirac points does not change with 



changing periodic field from B{r) to —B{r). Indeed, any 
gradual variation of the field amplitude from Bq to —Bq 
does not change the Chern numbers because, as one can 
see from the band structure calculations, this does not 
produce any field-induced band crossings. This does not 
depend on the value and sign of A and therefore is also 
valid in the limit of A — > 0. Since the sum of Chern 
numbers and the resulting AHE at /x = is invariant 
with respect to field inversion, then, due to the antisym- 
metry of AHE to the field, it should be zero. Thus, in 
the vacuum state with fj, = there is no Hall effect in 
graphene. 

If the chemical potential /i 7^ and is located some- 
where within the gap, the AHE is not necessary zero and 
can be calculated from the following formula, which gives 
the quantized values of the Hall conductivity 




Here Ch„ = ^ J d^k Jlkn is the Chern number of the 
n-th energy band, and the first sum in (14) runs over the 
occupied energy bands. 

Using the obtained results for the Chern numbers we 
find that the contribution from the /C point changes (in 
units of e^/h) from —1 to when the parameter A 
changes from negative to positive values. Similarly, the 
contribution from the /C' switches from to —1 when A 
changes from A < to A > 0. The sum of contributions 
from both K, and K,' is equal to —1 for any sign of A, and 
the sign of A determines, which of the points JC or JC 
is contributing to the AHE. But in the point of A = 0, 
which corresponds to the gapless model of graphene, the 
contributions from /C and JC' are exactly the same. There- 
fore, we come to conclusion that at A = each of these 
points gives the fractional value of quantized Hall con- 
ductivity equal to —1/2. 

VI. GAPLESS DIRAC MODEL 

Let us consider now the Hamiltonian (1) with A ~ 
0. It has the symmetry t^^Htz = — "H, which means 
that if Tpn is the eigenfunction of % with energy e„ then 
TzTpn corresponds to the state with energy — £„. The 
symmetry also means that the Berry phase calculated 
for any positive (-f n)-band is exactly the same as for the 
{—n) band. 

For any k-point including F we use the unperturbed 
wave functions (3). One can calculate the Berry phase 
7c along the contour of Fig. lb, excluding the F point 
and using 

7c = arg (r|M) {M\K) {K\T') (F'|r) , (15) 

where F and F' correspond to the states with k ~ 
along the T-M and F-A' lines, respectively. It gives us 
for the /C point and e > 0, 7c = — 7r/12. The contour 
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FIG. 2: The electron band structure of graphene (A = 0) for 
different intensities of the field: (a) a = 0.05, (b) a = 0.2, (c) 
a = 0.5, (d) a = f.O. 



around the Brillouin zone gives us 7+1 = 127c + 70, 
where 70 is the Berry phase for a small contour around 
the r point. The calculation gives 70 = tt. As a result we 
find the integer Chern number Ch+i = 0. However, this 
Chern number docs not determine the Hall conductivity 
corresponding to Eq. (14) because the contribution of 
7o should be excluded from consideration as it does not 
depend on the external field. Thus, to calculate the AHE 
we should use the Berry phase without enclosing F point, 
and we finally find that the contribution to AHE from 
any of the /C or /C points is equal to —1/2 in accordance 
with what we found before. 

Chern numbers 
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FIG. 4: Hall conductivity per one spin and one valley as a 
function of magnetic field for different values of electron den- 
sity: (a) no — 4, (b) no = 8, (c) no — 12, (d) no = 16, 
where no is the total number of electrons per elementary cell 
of graphene in periodic field. 



VII. NUMERICAL CALCULATIONS 

Using Hamiltonian (1) we performed numerical calcu- 
lations of the energy band structure for different intensity 
of the periodic field. In these calculations we used the 
real numerically calculated distribution of magnetic field 
created by the lattice of magnetic nanorods, not restrict- 
ing ourselves to harmonic approximation. The results are 
presented in Fig. 2 for different values of the dimension- 
less parameter of intensity of periodic field a, defined as 
a = Bqq'^ /4Tr4>o, where ^0 is the flux quantum. Figure 2 
demonstrates that the energy gap between the (-1-1) and 
(-1-2) bands appears at the field corresponding to a > 0.2. 

We also calculated numerically the Chern numbers as 
a function of the field intensity. They are presented in 
Fig. 3 for the first three bands. It turns out that the 
results of perturbation theory for the low-energy bands 
n = ±1 remains valid with increasing a. Figure 4 
presents the Hall conductivity per one spin and one valley 
as a function of concentration of electrons. The dashed 
line corresponds to the unfilled energy band, where an 
additional contribution from the Fermi surface should be 
taken into accouni^H but we do not consider this case. 



VIII. ESTIMATION OF PARAMETERS 



FIG. 3: Numerically calculated Chern numbers of the first 
positive energy bands for small A 7^ as a function of a: (a) 
Ch+i for A > 0, (b) Ch+i for A < 0, (c) Ch+2, (d) Ch+3. 



The results of numerical calculations show that the 
first gap between the (-1-1) and (-1-2) bands appears at 
a > 0.2. In the case of iron nanorod lattice with a lat- 
tice parameter a ~ 150 nm, we find that the magnitude 



of Bq at a distance about 15 nm from the nanorod lat- 
tice is Bq ~ 0.2 T— and that a ~ 0.2, which is in good 
agreement with the criteria of gap formation. 

One can also estimate the electron density n corre- 
sponding to the insulating state when the chemical po- 
tential is located in the first gap. This state is char- 
acterized by one electron per elementary cell, n = 4/iSo, 
where Sq = 0^^/3/2 and factor 4 takes into account de- 
generacy with respect to spin and valley in graphene. Us- 
ing the relation n — fi^ /ttv'^, where v/h ~ 1.5 x 10^ m/s 
is the electron velocity, we find ^ ~ 40 meV. All neces- 
sary parameters are quite reachable for the technology of 
graphene. 

One can estimate the characteristic energy Eq = 
hv/a ~ 40 meV for a = 150 nm. From the band struc- 
ture calculations we expect the gap to be of the order 
of Eq/G. So, the gap is ~ 90 K. For the Zeeman split- 
ting in graphene we have gfisBo — 0.5 K, so we expect 
that Zeeman splitting is negligibly small for this choice 
of parameters. 



IX. DISCUSSION 

We found that the AHE related to a single Dirac point 
in graphene is fractionally quantized, corresponding to 
quantum number —1/2. The quantization of AHE takes 
place if the chemical potential is located in the gap be- 
tween the bands (-1-1) and (-f2), which can be easily 
achieved with the existing experimental technique. The 
fractional quantization is related to the Dirac points, 
which is the monopole (for the K- point) or antimonopole 
(for JC') in the k space, but the Berry phase related 
to these points does not contribute to the AHE. Cor- 
respondingly, the measured AHE in graphene with two 
nonequivalent Dirac points and two spin orientations 
is integer for the same parameters, <7h — ne^/h with 
n ^ -2. 

Our consideration of the AHE is essentially based on 
the statement that in the clean case with /i = the AHE 
is zero. In the above-studied model with A 7^ 0, it is 
related to the invariance of Hamiltonian and Chern num- 
bers of the filled bands with respect to the inversion of 
field, which contradicts to antisymmetry of AHE. It does 
not depend on sign of A and leads to the zero vacuum 
AHE in the limiting case of A = in graphene. 

In this connection, it is important to compare our 
results with the result of Haldane's modeli^, in which 
the AHE can be nonzero even in the vacuum state. In 
the Haldane's model the periodic field has the same pe- 
riodicity as the honeycomb lattice, and the model in- 
cludes additional next-neighbor hoppings. Then the 
next-neighbor hopping integral in such a fast alternat- 
ing field acquires nonzero phase factor, which leads to 
the brake of electron-hole symmetry. Obviously, this 
model is not invariant with respect to the field inver- 
sion, B{r) —B{r). As a result, our arguments about 
vanishing of AHE in the ground state of graphene are not 
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FIG. 5: Brillouin zone and Berry sphere of the Haldane's 
model in tlie case of A < 3V3t2 sm 0. 



applicable to the model of Ref. [1^1 • 

To clarify this point we reconsider the Haldane's model 
by calculating the Hall conductivity as the Berry phase 
of electrons in the lower subband along the path encom- 
passing the Brillouin zone of the honeycomb lattice. The 
Hamiltonian of Ref. [l5j can be presented as 



Hh = e(k) + Act • n(k) 



(16) 



where the unit vector n(k) parametrizes 2x2 pseudospin 
Hamiltonian. Here we denoted 

e(k) = 2t2cos0^cos(k • bj), (17) 

i 

n,(k) = (ii/A)^cos(k-a,), (18) 

i 

nj,(k) = (ti/A)^sin(k-a,), (19) 

i 

n^{k) = [A-t2sin(?!)^sin(k-bj)]/A, (20) 

where the lattice vectors a,; and are defined in [Tsj 
and the parameter A is introduced to provide |n|=l. The 
parameter ^2 is the next-neighbor hopping integral and 
the phase is due to this hopping in the field with the 
periodicity of crystal lattice. 

The Berry phase for the motion along the path 
(0, ki, kj\/, — ka, 0) is equal to the spherical angle of the 
mapping k — > n(k) to the Berry sphere (see Fig. 5). Us- 
ing (17)-(20) we find that in the case of A > 3\/3t2 sin0 
the vector n(— ka) points down coinciding with n(ki), 
which leads to the Berry phase 7c = and Chern num- 
ber Ch = 0. In the opposite case of A < 3^/3t2Sm(f>, 
vector n(— ks) points up (as shown in Fig. 5) lead- 
ing to = — 7r/3 and, correspondingly, we obtain in 
this case Ch — —1. The same result follows from the 
explicit calculation of the Berry phase using relation 
7c = arg (0|ki}_(ki|kM) (kA/l - kg) (-kajO). This is the 
result of Ref. 115] found by geometrical method. 

When the constant of field lattice a is much larger than 
the lattice constant of graphene oq, the electron- hole and 
inversion symmery of graphene should be restored. In 
this limit, the field flux through an elementary cell is 
nonzero breaking locally the inversion symmetry but the 
number of cells with positive and negative fiux is equal 
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at the scale I ^ a. Besides, in any realistic structure, the 
spatial fluctuations of field periodicity are much larger 
than the lattice constant, which makes it impossible to 
provide the compatibility of two different lattices. 

We believe that the quantized Hall effect in the struc- 
ture with graphenc! on top of the magnetic nanolattice 
can have important practical applications similar to the 
usual quantum Hall effect in magnetic field. 
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